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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of GSC 6214- 
210 A and B, a solar-mass member of the 5-10 Myr Upper Scorpius association with a 15 ± 2 Mj up 
companion orbiting at ~330 AU (2"2). Previous photometry and spectroscopy spanning 0.3-5 /im 
revealed optical and thermal excess as well as strong Ha and Pa /3 emission originating from a circum- 
substellar accretion disk around GSC 6214-210 B, making it the lowest mass companion with un¬ 
ambiguous evidence of a subdisk. Despite ALMA’s unprecedented sensitivity and angular resolution, 
neither component was detected in our 880 /im (341 GHz) continuum observations down to a 3-<r limit 
of 0.22 mJy/beam. The corresponding constraints on the dust mass and total mass are <0.15 
and <0.05 Mj up , respectively, or <0.003% and <0.3% of the mass of GSC 6214-210 B itself assuming 
a 100:1 gas-to-dust ratio and characteristic dust temperature of 10-20 K. If the host star possesses a 
putative circum-stellar disk then at most it is a meager 0.0015% of the primary mass, implying that 
giant planet formation has certainly ceased in this system. Considering these limits and its current 
accretion rate, GSC 6214-210 B appears to be at the end stages of assembly and is not expected to 
gain any appreciable mass over the next few Myr. 

Subject headings: accretion disks — brown dwarfs — stars:individual (GSC 6214-210) 


1. INTRODUCTION 

Giant planets are common products of protoplane¬ 
tary disk evolution but many of the details of their 
formation remain obscure and untested. Core accre¬ 
tion is the dominant formation channel at small sep- 


arations of <10 AU dPollack et al 


12005 b while disk instability dBossI 


199 


lAlibert et all 

and turbulent 

fragmentation (fBatei (20091) operate on wider orbits at 
tens to hundreds of AU. In these outer regions (~10- 
300 AU), high-contrast direct imaging surveys are re¬ 
vealing that the tail end of the substell ar com panion 
mass function is likely to be continuous ■: Brandi ii al 
2014) , extends down to at l east «5 Mj up (jRameau et al 


^ U.U VV 11 \JVJ OpU IGUiOU ' —' ± VJL J pi p 

20131 iKuzuhara e~a l. 2013), and ove rlaps with the most 


massive planets lik e HR 8799 bc de dM arois et all I 
iMarois et al.l l2010h and /3 Pic b ((Lagrange et al .112001 
which probably formed in massive protoplanetary disks. 

In all formation scenarios, giant planets accrete most 
of their mass from circum-planetary disks, which natu¬ 
rally form around protoplanets in an analo gous way to 
circumstellar disks around protostars (|Vorobvov fc Basul 
120101 1 Ward fe Canupll2C)10D . Circum-planetary disks reg¬ 
ulate angular momentum transport and accretion as pro¬ 
toplanets grow, u ltimately determining the final masses 
of gi a nt planets dQuillen fc Trilling! 119981 iLubow et all 
119991 iCanuD fc Wardl 120021 iZlml 120151) . At the earli- 
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est stages of giant planet assembly (<3 Myr), accretion 
shocks from circum-planetary disks onto protoplanets 
help to define planets’ initial entropy levels and which 
evolutionary pathwa ys ( hot, warm , or cold start mod¬ 
els) they will follow ilMarlev et al.l[2007h . Moreover, the 
formation, bulk composition, and orbits of moons are 
also g overned by properties of circum-planetary disks 
dHeller et al.ll20l4) . Physically, these disks are expected 
to be quite thick ( h/r > 0.2), possess little flaring com¬ 
pared to their circum-stellar analogs, and be tidally trun¬ 
cated to ~l/3 i?Hin from gravi t ational influence of the 
host star (lAvliffe fc Ba.tel 120091 iMartin fe Lubowl 120111 : 
iShabram fc Bolevll2013lh 

This epoch of planet growth is challenging to directly 
study because of the small angular resolution and high 
contrasts required. The possible discovery of extended 
thermal exce ss emission from the ca ndidate protoplanets 
LkCa 15 b (1 Kraus fc Irelandl 120121) and HD 100546 b 
(iQuanz et al’l I2013T ) may represent the first direct de¬ 
tections of circum-planetary disks, though the masses 
and immediate environments of the se embedded com¬ 
panions are still poorly constrained (jlsella et al.l [2014 1. 
The remnants of circum-planetary disks are evident from 
the prograde, regular satellites orbiting gas giants in 
our Solar System as well as the (possibly transient) 
ring systems surrounding Saturn, the subst ell ar com- 
pankm 1SWASP J140747.93 -394542.6 dMamaiek et all 
l2012|: pCenworth y et al.ll2014fl . and perhaps Fomalhaut b 
(iKalas et al.ll200fih 

Over the past decade a population of young 
(<10 Myr) planetary-mass companions on ultra-wide or¬ 
bits (>100 A U) has been u ncove red with direct imaging 
(see Table 1 of lBowler et al.ll20T4i ). Many of these possess 
their own accretion disks, offering unique opportunities 
to study the structure, div ersity , and evolution of circum- 
planetary disks. Recently iKraus et al.l (120151) presented 
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Fig. 1.— Spectral energy distributions of GSC 6214-210 A (red) and B (green). GSC 6214-210 B exhibits optical and thermal excess 
along with extraordinarily strong Ha; emission. The unresolved pair shows slight excess beyond 12 ^m that could originate from the 
circum-substellar disk around GSC 6214-210 B and/or a low mass disk around the primary. Our ALMA 880 fn n 3 -a upper limit (blue) 
rules out substantial disks from both components. Photometr y for the primary are from CMC-14 (r'), DENIS ( UK ), 2MASS ( JHK ), 
WISE (3.4, 4.6, 12, a nd 22 ^m), and Sp itzer /MIPS (24 um; IBailev et all 120131) . Note that the 3.1 /im, 3.3 pm, and M' photometry of 
GSC 6214-210 A from [Bailey et ahl (120131 ) are inferred rather than independently measured. Photometry for GSC 6214-210 B are from 
llreland et all BOTH. JHKL'j,*S ailev et aH (1201 31. 3.1 um , 3.3 pm, M'), and IZhou et al.l (120141 , 0.33-0.85 pm). The BT-Settl 4200 K/4.0 dex 
and 2400 K/4.5 dex models (dashed lines; [Allard et al. 2010]) scaled to the J-band photometry are shown for comparison. 
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1.3 mm Atacama Large Millimeter/submillimeter Array 
(ALMA) continuum observations of FW Tau C, a young 
(«2 Myr), wide-separation (330 AU) companion with a 
possible mass as low as 6-14 Mj up dKraus et all 120141) . 
The clear detection at mm wavelengths implies a dust 
mass of 1-2 Mg, sufficient to form a system of rocky 
moons. However, the high veiling and possible edge-on 
architecture of the this accretion disk makes the mass of 
FW Tau C highly u ncerta in an d poss ibly well above the 
planetary regime dBowler et al.ll20l4l) . 

Here we present ALMA 880 /im continuum observa¬ 
tions of GSC 6214-210 B, a «15 Mj up companion or¬ 
biting the young (5-10 Myr) Sun-like star GSC 6214- 
210 A in Upper Scorpius (USco) at 330 AU. With a mass 
at the deuterium-burning limit and abundant evidence 
of ongoing accretion, GSC 6214-210 B is the lowest- 
mass companion with a circum-substellar disk currently 
known. ALMA’s unprecedented sensitivity and angu¬ 
lar resolution at sub-mm and mm wavelengths presents 
a unique opportunity to characterize a sub-disk at the 
boundary between the brown dwarf and planetary mass 
regimes. Additionally, ALMA observations provide im¬ 
portant contextual information about any circum-stellar 
disk surrounding GSC 6214-210 A, or lack thereof in case 
of a non-detection. 

2. THE GSC 6214-210 AB SYSTEM 

GSC 6214-210 is a well-established member of the USco 
star-forming region (IPrei bisch et al.l 1 19981 llreland et all 
120111: ILuhman fe; Mamaiekl 120121) . Low-mass evolution¬ 
ary models point to an age of ~5 Myr for this complex, 
though more rec ent tu r n-off ages su ggest it is somewhat 
older (sslO Mvr: iPecaut et aLll2012D . Based on its spec¬ 
tral type (I<5 ± 1) and modest reddening (Ay~0.6 mag), 


evolution ary models imply a mass of 0.9 ± 0.1 Mp, for the 
primary (IBaraffe et al.l 119981 : IBowler et al.l 1201 ill . GSC 
6214-210 is an evolved Class III object with no signs of 
a substantial disk out to 12 /im. 

The low-mass companion to GSC 6214-210 was first 
i dentified at 2"2 (330 AU) from Keck/AO imaging by 
IKraus et all (12008) and confirmed to be comoving with 
the primary bv llreland et all (|2011l) . Evolutionary mod¬ 
els imply a mass of 15 ± 2 M Jup from its age (5- 
10 Myr) and l umino s ity (lo g Lboi/R© = -3.1 ± 0.1 dex). 
IBowler et al.l (| 2011 l 120141 ) found a near-infrared spec¬ 
tral type of M9.5 ± 1 for GSC 6214-210 B as well as 


strong Pa/3 emission at 1.282 /im (EW=-11.2 A), indi¬ 
cating vigorous accretion from a circu m-subste l lar dis k. 
Recen tly, 1.0-2.4 /im spectroscopy bv lLachapelle et all 
(120151) also revealed Pa/3 emission as well as weak Br 7 
emissio n at 2.166 /tm. Optical imaging with HST by 
I Zhou et~aLl (2014) uncovered extraordinarily strong Ha 
emission (EW ~-1600 A) and continuum excess short- 
ward of ~6000 A, corresponding to an accretion rate of 
M~1Q ~ 10 - 8 Mp) yr~h 3-5 /tin imaging bv llreland et all 
(I 2011 D and IBailev et al.l (I2013D supports the presence of 
a warm disk from thermal excess emission above the 
photospheric level. Unresolved 22 /im photometry of 
GSC 6214-210 AB from WISE and 24 /im photometry 
from Spitzer ([Bailev et al.l 120131 1 reveals a slight (but 
significant) excess, which could be from an evolved disk 
around the primary and/or a massive disk surrounding 
the companion (Figure |Tj). Synthesized IU4-band pho¬ 
tometry of the 4200 K/4.0 dex and 2400 K/4.5 dex BT- 
Settl models in Figure [T| implies a range of flux densities 
for the companion between 9.92xlO -20 W m ~ 2 proT 1 
(or 2.19xl0” 18 W m -2 at 22 /im, assuming no excess) 
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Fig. 2.— ALMA 880 fim continuum image centered on GSC 6214- 
210 AB. The inset shows a 2 // x4 // A'-band Keck/N IRC2 adaptive 
optics image of the pair from [Ireland et al.l ^20111 ). The same re¬ 
gion is depicted as a dashed box in the ALMA image. Neither 
component is detected with a 3-cr upper limit of 0.22 mJy (assum¬ 
ing a point source). The field of view is 10"xl0" and the 0 / /3x0 / /4 
beam FWHM is shown as a white ellipse in the bottom left corner. 
The image has been corrected for the nonuniform sensitivity of the 
primary beam across the field of view, and the color stretch spans 
=b5cr (d=0.37 mjy/beam). 

to 1.2 ± 0.6xl0 -17 W m- 2 /mr” 1 (or 2.6 ± 1.3xl(r 16 
W m -2 at 22 /im, assuming the excess is entirely at¬ 
tributed to the companion). 

3. ALMA OBSERVATIONS 

GSC 6214-210 AB (a=16:21:54.67, £=-20:43:09.2) was 
observed with ALMA in Band 7 at 880 /im (341 GHz) 
during Cycle 2 Early Science operations on UT 2014 June 
30. 31 12-meter antennas targeted this system with the 
dual-polarization setup in four spectral windows centered 
at 334.0, 336.0, 346.0, and 347.8 GHz. The 346.0 GHz 
window encompases the CO (J=3-2) line and was sam¬ 
pled in 1920 channels each with a width of 0.98 MHz (fre¬ 
quency division correlator mode), while the other three 
continuum windows comprised 128 channels at 15.6 MHz 
each (time division correlator mode). The total on-source 
integration time was 11.6 minutes and the total effective 
bandwidth was 7.7 GHz. The maximum array baseline 
reached 650 m, resulting in a synthesized beam FWHM 
of 0 / .'32x0 , .'41 at a position angle of 96.3° and a field of 
view spanning 18"xl8". 

The visibility data reduction and image reconstruc¬ 
tion were carried out with CASA version 4.2.2. The 
primary calibrator was the nearby quasar QSO J1625- 
2527, which was targeted intermittently during the sci¬ 
ence observations. No line emission was detected in the 
346.0 GHz data so this window was incorporated with 
the continuum regions to create an averaged continuum 
image. No detections are present above the background 
rms level of 0.074 mJy/beam as measured in the central 
2!!2x2!!2 region. Finally, the continuum image was cor¬ 
rected for the non-uniform sensitivity across the field of 
view by dividing by the primary beam response. The 
result is shown in Figure [2] 


4. RESULTS 

There is no significant emission above the background 
level (0.074 mjy/beam rms) at the location of GSC 6214- 
210 A or its companion, implying a 3 -a upper limit of 
0.22 mJy (assuming a point source). With a spatial 
resolution of ss0'/4, both components would be easily 
separated in our observations. Assuming the contin¬ 
uum emission is optically thin, an upper limit on the 
dust mass (M<j ust ) can be derived from our ALMA flux 
density constraint, the distance to the system (145 pc), 
a characteristic isothermal dis k tempera ture (Tb ), and 
a dust opacity («„) following iHildebrandl (11983f l. For 
consistency with previous work we use the frequency- 
dependent dust opacity relationship k v = 10(i//10 12 Hz) 
cm 2 g” 1 from iBeckwith et al.l (119901) . At 880 /im, the 
dust opacity Kssoum is 3.4 c m 2 g~ . 

Following I Andrews et al.l (120131 ). we adopt the dust 
temperature-stellar luminosity relationship Tc = 25 
(. L/Lq) 1//4 for the primary GSC 6214-210 A. This yields 
Tc =20 K using the lumi nosi ty of log L/Lb 0 \ = ~ 
0.42 ± 0.08 dex from iBowler et al.l (j2011[ ). Heating 
of GSC 6214-210 B’s outer disk is probably domi¬ 
nated by GSC 6214-210 A or diffuse interstellar radia¬ 
tion rather than irradiation by GSC 6214-210 B itself. 
We therefore assume a dust temperature between 10- 
20 K for GSC 6214-210 B, which is within the range of 
disk mid plan e te mperatures for isolated brown dwarfs 
(|Ricci et al.ll2014fl . 

The 3-cr dust mass upper limit is <0.05 M® (3.9 lunar 
masses) for the primary and <0.05-0.15 M® (3.9-12.6 
lunar masses) for the companion. For a gas-to-dust ra¬ 
tio of 100, this implies putative disk masses of <5 M® 
for GSC 6214-210 A and <15 M® for B. These results 
are most sensitive to the characteristic disk tempera¬ 
ture, dipping to 0.03 M® (2.8 lunar masses) of dust for 
a slightly warmer temperature of 25 K, and increasing 
to 0.9 M® (76 lunar masses) for a temperature of 5 K. 
The lower limit on dust temperature is set by the diffuse 
interstellar radiation field and cosmic rays, which even 
for dense starless cloud cores that are shielded from ex¬ 
ternal radiation is only about 10 K. So although our dust 
mass upper limit is highly sensitive to the dust temper¬ 
ature, we consider the dust mass of ~0.15 M® set by a 
Tc =10 K to be a fairly strict upper limit on sub-rnm- 
sized particles and smaller. The total disk mass also de¬ 
pends on the grain size distribution, which is especially 
important for an older, presumably evolved system like 
this. If planetesimals have formed then they can also lock 
up a significant amount of solid mass, though there is an 
upper limit to this storage since fragmenting collisions 
would grind them back down to ~mm sizes that would 
be detectable with these ALMA observations. Moreover, 
since the gas-to-dust ratio in protoplanetary disks is gen¬ 
erally highly uncertain and evolves over time, the total 
disk masses in this work are approximate at best. 

5. DISCUSSION 

Our deep 880 /tm upper limit implies that GSC 6214- 
210 B’s disk has a low mass (<0.05 Mj up ) and lacks 
a significant amount of cold dust in its outer regions. 
This is broadly consistent with expectations from circum- 
planetary disk simulations, which find that subdisks are 
truncated at ~l/3 Rum as a result of the gravitational 
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Fig. 3. — Top: Sub-mm/mm dust mass measurements as a func¬ 
tion of stellar mass for the Myr Taurus (blue) and «5—10 Myr 
USco (orange) star-forming regions. Triangles denote 3-<r upper 
limits. Our ALMA constraints for GSC 6214-210 A and B are 
over an order of magnitude lower than f or most previous sur ¬ 
vey s in these regions. Tauru s data are from [Andrews ct al. ( 2013) 
and [A kcson & Jensen ( 2 014]); USc o data are from Mathew s et al.l 
<20lf) and [Carpenter ct al. (20f4|). Dust masses from Andrews et 
al. are derived from their 890 fi m flux density measurements and 
dust temperature-stellar luminosity scaling with a lower boundary 
of 10 K adopted for the lowest-luminosity brown dwarfs in their 
sample. Ordinates on the right-hand axis denote total disk masses 
assuming an interstellar-like 100:1 gas-to-dust ratio. For compar¬ 
ison, t he dust mass measurement of FW Tau C from IKraus et al.l 
(120151) is labeled. Bottom: Similar to the upper panel, but show¬ 
ing dust mass to stellar mass ratio (and total disk mass to stellar 
mass ratio) as a function of stellar mass. In terms of mass ra¬ 
tio, our non-detection of GSC 6214-210 B’s disk is not particularly 
unusual. 

influence of the host star dShabram fe Bolevll2013lh For 
GSC 6214-210 B, this corresponds to a radius of ~19 AU. 
Interestingly, the free-floating brown dwarf OTS 44 in 
Chamaeleon I has a similar mass (rU 2 Mj up ) and accre¬ 
tion rate (~8xl0 -12 M 0 yr _1 ) to GSC 6214-210 B, but 
possesses a significantly more massive disk of ~30 M 0 
and is co nsistent w ith h aving a spatial extent out to 
~100 AU (jJoergens et al.1 2013 1. These differences could 
be caused by disk evolution owing to the different ages 
of these objects (~2 Myr versus 5- 10 Myr), environ¬ 
mental influences from the nearby host star GSC 6214- 
210 A, natural variations in the physical properties of 
circum-substellar disks, and/or differing formation mech¬ 
anisms. It is also possible, though more speculative, 


that GSC6214-210 B originally formed closer in and 
was dynamically scattered to a wide separation from an 
as-yet-undiscovered more massive companion. In this 
scenario, GSC6214-210 B’s disk could have been trun¬ 
cated, explaining its hot, inner component but low over¬ 
all mass. We note that it is possible the disk may be very 
small (<0.4 AU) and optically thick given the bright¬ 
ness temperature limit of <0.02 K, but that would re¬ 
quire long-term stability of the disk in the face of short 
viscous timescales. Ultimately, the growing population 
of plan etary-mass compan ions being found at hundreds 
of AU (|Kraus et al.l 120ill combined with ALMA’s un¬ 
precedented sensitivity and spatial resolution will make 
it possible to explore these effects for larger samples of 
free-floating and bound planetary-mass objects spanning 
ages of 1 10 Myr. 

In Figure [5] we compare our ALMA constraints to pre¬ 
vious sub-mm/mm surveys in the Taurus and USco star¬ 
forming regions. Our upper limit for GSC 6214-210 A 
indicates there is a wide range of outcomes for proto¬ 
planetary di sks around Sun - like st ars i n USco. Larger 
progra ms bv IMathews et al.l (|2012f> and ICaroenter et all 
(2014) found disk masses in excess of 10 Mj up , implying 
a spread of over 3 orders of magnitude in disk mass when 
taking into account our non-detection. Indeed, if the pri¬ 
mary star GSC 6214-210 A harbors a disk, it must be at 
most a putative 0.0015% of the host star mass— much 
lower than the median dis k-to-star mass ratio o f ~0.3% 
found at younger ages by (Andrews et al.l (1201311 . Alto¬ 
gether it is clear that giant planet formation has ceased 
in this system. 

Similarly, there is not enough dust in GSC 6214-210 B’s 
disk for gas giant “moons” to form, but we cannot 
rule out future rocky moons with mas s es of »0.15 M 0 
(«12 lunar masses). ICanup fe Wardl (12006H find that 
the shared mass ratio of Jupiter, Saturn, and Uranus 
to the total mass of their regular moons (~10 4 ) can be 
explained by a balance between inflowing material into 
their circum-planetary disks, which provides the raw ma¬ 
terial for creating new moons, and satellite loss through 
orbit decay. In this scenario the current regular moons 
around these gas giants are simply the last generation 
of those formed from their subdisks. The current dust 
content of GSC 6214-210 B is a factor of at least 3 lower 
than this planet-moon mass ratio, so if this phenomenon 
is universal then any moons already formed or in the pro¬ 
cess of forming in this system will probably be the last 
ones. 

Assuming a steady-state accretion rate, the lifetime 
of GSC 6214-210 B’s disk can b e estimate d from its 
accretion rate (10 -10 ' 8 M 0 vr -1 : [Zhou et al.l I2014T) and 
disk mass constraint. At most, accretion will continue 
for another ^3 Myr, though the fractional gain in planet 
mass is only <0.4%. The GSC 6214-210 AB system 
therefore appears to be in its last stages of assembly, 
with giant planet formation having ceased around the 
primary star and moon formation probably in its final 
stages around the companion. 


We are grateful to the referee for helpful comments, 
Jonathan Swift for productive discussions about pur¬ 
suing this idea, and Vanessa Bailey for providing 
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zero point flux densities for MMT and LBT filters. 
This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2013.1.00487.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRC (Canada) 
and NSC and ASIAA (Taiwan), in cooperation with the 
Republic of Chile. The Joint ALMA Observatory is oper¬ 
ated by ESO, AUI/NRAO and NAOJ. The National Ra¬ 
dio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. We utilized data 


products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronau¬ 
tics and Space Administration and the National Science 
Foundation. NASA’s Astrophysics Data System Biblio¬ 
graphic Services together with the VizieR catalogue ac¬ 
cess tool and SIMBAD database operated at CDS, Stras¬ 
bourg, France, were invaluable resources for this work. 

Facilities: ALMA 
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